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УДК 681.51 
О.И. АТРОЩЕНКО 


СИНЕРГЕТИЧЕСКИЙ СИНТЕЗ УПРАВЛЕНИЙ 
ДЛЯ НЕЛИНЕЙНОГО ОБЪЕКТА УПРАВЛЕНИЯ 


Рассматривается применение метода АКАР (Аналитическое конструирование аг- 
регированных регуляторов по заданным инвариантным многообразиям) синергети- 
ческой теории управления для синтеза законов управления нелинейными распре- 
деленными объектами управления. Объектом управления выбран водоносный гори- 
зонт центрального участка Кисловодского месторождения минеральных вод. 
Ключевые слова: коэффициент фильтрации, граничные условия, математическая 
модель, водоносный горизонт, фазовое пространство, инвариантное многообразие, 
переходной процесс, управляющее воздействие, фазовый процесс. 


Введение. Водоносный горизонт — это объект с распределенными пара- 
метрами, поведение такого объекта существенно зависит не только от вре- 
мени, но и от пространственных координат [2, 4]. На центральном участке 
Кисловодского месторождения расположены пять эксплуатационных и пять 
контрольных скважин; эксплуатируется Валанжинский ВОДОНОСНЫЙ гори- 
зонт [2]. 
Постановка задачи. Имеется математическая модель системы управле- 
ния, необходимо определить законы управления методом АКАР, обеспечи- 
вающим минимальное истощение эксплуатационных запасов подземных 
вод, т.е. снижение уровня водоносного горизонта [2]. 
Методы испытаний. Нелинейная математическая модель системы управ- 
ления в переменных состояния имеет вид [2]: 

д, = - ах + Ю4х, - Вах, + Кох - Кзх, + ких, - (Е - зи (х,)) 
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х, = ках, - 2юах, + Ка х, - Ка х, + КЧх. - Ках, + КоЧхь- 
4%. 
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ж= Юах.- 2 аж + Кам + Юж - 2х - Ка ж + Каха - 


- (Е - 15))—^®; 
Ах 


х = аж - ах + КоЧ хи + ах. - 2х - Ко. хо + Кох - 





- (Е - 15) ; 
Ах 


Хо = Юм - Ко хо + ах, - 2 о хь- Кох + Кро Хоо ; (1) 
Хи = - Ка + Кам» - К хи+ Ко хь + Вах, - КФХ - 
- (Е - это) 7; 
Ах 


Хр = Ка хи - 2 ахь + Кзахь - Кох + КФ хи + ках, - КоЧхь- 


- (В, - чт, уе; 
Ах 


Хз = А хь - 2 зах + Каха - КЗЧ хь + Юз + ЮЧух, - Кз@хь- 
-(В,- зто); 
Ах 
Ха = Аза хз - 2. ха + КФХ - Каха + Кох + К азх, - Каха - 


- (Е - зто) в, 
Ах 


Хх; = Каха - Кз@хь- КБ Кох» + Юж - 5х5 ; 
Хб = - Юва + зах + К хи- 2+ Юоужь - КеЧхь- 
- (В - тб) и; 
Ах 
Хи = - Аб@хь - Каха + Юз ха + Кох - 2х + КЧух, - Кох - 
-(Е,- вт.) 28 _ т: 
Ах 
Хз = Кама - 2Азахз + Коахо + Юз хз - Аза + уж - Юз@зхи - 
Е вш(х,)) 9—7 и, ; 
Ах 
Хо = АвЧхь- 2КоЧ хо + Кох + Каха - Кох + Кох - Кох - 
- (Е - ть) и, ; 
Ах 


Хр = Комо - Кох + 54 м5- Кох» + Кох - Ко Ч5Хьо - 5. 

Система (1) является многомерной, нелинейной, многосвязной про- 
странственно-неинвариантной моделью, имеющей пять каналов управле- 
ния. Коэффициенты фильтрации А - А», зависят от пространственных коор- 
динат, скорость движения водоносного горизонта Ё= К, - т(х,) является 
нелинейной величиной, А = 1.5 м/сут. 

В систему (1) входят следующие параметры и коэффициенты: 
х -х» - переменные состояния, физический смысл которых — понижение 
уровня водоносного горизонта в точках дискретизации распределенного 
объекта управления. Состояния х- х; определяют понижение уровня во- 
доносного горизонта в контрольных скважинах, а состояния х5- х» озна- 
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чают понижение уровня водосного горизонта в точках расположения экс- 
плуатационных скважин; и,- и; — водозабор (дебит) из эксплуатационных 
скважин, физический смысл которого — управляющее воздействие на 


объект управления; ! — коэффициент упругоемкости коллектора в 1/м. 
и Но Ло 
Е а 2 — шаги дискретизации по осям 
Хх, у, 2. 

Согласно методу АКАР, необходимо синтезировать законы 
управления и, и,, и., и., и; как функцию координат состояния объекта 
и(х), обеспечивающие минимальное понижение уровня в контрольных 
скважинах на один метр [2, 5]: 

О еее ва, ее, Ее. 22) 

В соответствии с числом независимых каналов управления может 
быть назначено пять инвариантов (стабилизация понижения уровня в точ- 
ках расположения контрольных скважин х,, х., х;, х., х.). Исходя из по- 


ставленной задачи (2), целесообразно выбрать параллельную совокупность 
внешних инвариантных многообразий: 





У „(№».--5)=0, т= 1,....5 
Ут = 6+9, ; = х+95; Уз= М1 0, 
У = №9194, И5- + Ф.. (3) 


Инвариантные многообразия (3) должны удовлетворять решению 
системы функциональных уравнений АКАР [1]: 


Ту уу, =0; ТУ @чу,=0; Ту. чу, =0; 


учу: =0; Ту (чу: =0. (4) 
При попадании изображающей точки системы на пересечение 
многообразий У,=0, у,=0, у,=0, у.=0, у.;=0, согласно принципу 
«расширения — сжатия фазового пространства», происходит декомпозиция 
исходной системы (1), и динамика системы описывается системой диффе- 
ренциальных уравнений пониженного порядка [1, 5]: 
х, = - ках, + ках», - ках, + 64. хе, - ках, + ах - 
- (В - ми, Убе, 
Ах 
х„ = ках, - 26 ах», + Вам, - Юаух,, + ках, - ках», + Ах, - 
-(В- щбь, ^^, 
Ах 
Ху = ках», - 25а, хз, + 4х - ах» + Ах, - ах, + зах, - 
- (Во - этбь, ут, 
Ах 
Хи Е ах, - 2х, + зах; - ах + 4х, - Ах, + Каха - 
- (В - зи, у, 
Ах 
Ху = ах, Е Ка, х., Е ху + Кох в ках. + 5х, 5, р 
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Хы = - Кахь + Ках, + Вах, - 2, хе, - КЧухь, - КьЯ9 1 - 


х,-х 
- (К - мии 
Ах 
ху = Ка - 2Кахь + Бажь, + КЧх, - 2х, - Юазху, - К ›- 
=: 
-(В,- зб, Ти, 
Ах 
Ху, — ках, - 2х, + ах, + 4х, - 2х, - АЧХ, - К 39.0. - 
х,-х 
- (у зи, )) м; 
Ах 
№, = аж, - 2. Ч, хо + оао + аж, - 2х, - Ч хо, - Ко 4 - 
хи -х 
- (Е, - эх» ))———*; 
Ах 
Хи 5 Кох, ы Кох + К х;, ы 2х, Е Кох, - Ко 36 5; (5) 


Ху а Ка х у + Ка, Хру < ах, - 264.9, + ках, ы К Чхиу и 


: Хр Ми. 
- (В - заб )) И, 
Ах 


Ху = ах - 2 ах, т ках, г К Ах, - К.4,9 ›+ ах. - Ка ху ы 


: Мзу Му , 
5 (В = (Ху )) " " ' 
Ах 


Ху = Ка, хь - 2 ЗЧ х,з + ах - Аз - Юз 4,9 3+ ах; - зах, > 


, М Му. 
- (В - (Аз —_—; 
Ах 


Ху = зах, - 24а ху + Аза х,5, - 4х, - Ко, 4+ Аха, 2 Каха, - 


: 5 ЛМ , 
- (о - зи, )) Е. 
Ах 


Ху = ках з К 5ах,5, - 2; 4>х,5, - 24.9 5+ ках, С ках, ; 
неизвестные функции #1, #5, 6%, 94, 9; - «внутренние» управления. 
На втором этапе синтеза для поиска внешних управлений ВВОДИТСЯ 


вторая группа «внутренних» инвариантных многообразий по числу каналов 
внутренних управлений: 


у „(№»:--№)=0, т = 6,...,10; 
6 -= №185; у = ху; Уз= №19: 
Уо= №+0о, Ую= Хю+ Фр, (6) 


где 9.-9, - неизвестные «внутренние» управления, которые необходимо 
доопределить в процессе синтеза. 
Инвариантные многообразия (6) должны удовлетворять решению 
системы функциональных уравнений АКАР: 


У @+У = 0; Ту (О+у,=0; 
Туту: 0; Ту о(+у,=0, ТУ 5 ()+Ую=0. (7) 
На пересечении совокупности внутренних инвариантных многооб- 
разий У‹=0, у,=0, у, =0, у.=0, У, =О динамика системы определяет- 
ся декомпозированной системой пониженного порядка: 
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х, = - ах, + Бах», - Юахь - 4,9 5- Казхь + Км - 


. х,-х 
- (Е - эта, )———*®; 
Ах 


Хх = ках, - 2, ах», + ках, - Ка. х», - Ка .- Юаух», + Кох г 


ж, -х 
- (В, - ше, )— (8) 
Ах 
ху = ЮЧх, - АЧХ, + Каха, - Кох, - №4. в - ЮзЧужу, + КзФуХ - 
хх -х 
- (Е, - пов, )— 5; 
Ах 
ху = Бах, - ах, + Бхь, - Юж, - Ко оо КаЧуха, + КФХ, - 
ху -х 
- (В - зто); 
Ах 


ху 2 К. в К, х., ы КЧ.х;, - 99 и- ках. + зах, ь 

На третьем этапе декомпозиции для поиска «внутренних» управле- 
ний вводится вторая группа «внутренних» инвариантных многообразий 
у, =0О, у, =0О, у: =0О, У .=0, у, =О по числу каналов с внутренними 
управлениями [1]: 

Ун=м-м=0; УЕ = О: Уз=%-%=0; 

Е Же 10 УЕ ж-%=0. (9) 

Инвариантные многообразия (9) должны удовлетворять решению 
системы функциональных уравнений АКАР: 

ТУ нчун=0; Ту р+ур=0; Туву =0; 

Ту а+У д = 0; Пу вуз =0. (10) 

В обратном порядке, совместно решая системы уравнений (8), (9), 

(10), определяем группу «внутренних» управлений Фз, 6%, 6%, во,вь. По- 
сле этого из совместного решения систем уравнений (5), (6), (7) определя- 
ем группу «внутренних» управлений @®,, ®., ©}, 9.4, $;, а из совместного 
решения систем уравнений (1), (3), (4) находим неизвестные внешние 
управления №, и. , из, и, и.. 
Результаты моделирования. Численное моделирование проводилось 
при следующих параметрах настройки регулятора: ФТ = 0.8; Т, = 0.9; 
Т.= 0.7; Т,=0.8; Т=0.9; Т,=0.8; Г=0.7; №=0.6; Т=0.9; Ть= 0.8, 
Т, = 0.9; 1, = 0.8; Т, = 0.9; Тл= 0.8; Т, = 0.8, а также следующих парамет- 
рах объекта управления: 1 = 0.000004 1/м; дх= 250; Ду-= 166; 42= 40; 
К = 0.05; №=0.06; №=0.06; №=0.07, №=0.07; №=0.08; К=0.02; 
Аз = 0.05; №= 0.01; Ж,= 0.01; №, = 0.01; &,= 0.12; й,= 0.02; К, = 0.09; 
Юз = 0.07; Ас = 0.07; К. = 0.02; Ку = 0.01; Ко= 0.11; А= 0.03 м/сут; В = 2.5 
. Начальные условия по переменным состояния были выбраны нулевыми: 
отсутствие внешнего воздействия. 

На рис. 1 представлены результаты моделирования: переходные 
процессы переменных состояний х,, х,, х;, х., х; по каждой контрольной 
скважине, стабилизирующихся к установившемуся значению (2). На рис.2. 
изображены изменения (стабилизация) во времени управляющих воздей- 
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ствий м, и,, и;, и., и; по каждой эксплуатационной скважине. На рис.3 
дана проекция фазового пространства на координаты х, - х, - х», фазовый 
портрет, пересечение внешних инвариантных многообразий У, = 0, у,=0 
в трехмерном фазовом пространстве, в котором фазовые траектории, стар- 
туя из произвольных начальных условий, под действием управлений попа- 
дают на пересечение многообразий, а затем движутся к установившемуся 
состоянию. На рис. 4 изображены переходные процессы дебитов эксплуа- 
тационных скважин в м? /сут. 
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Рис. 1. Переходный процесс Рис. 2. Управляющие воздействия 
переменных состояний 
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Рис. 3. Фазовое пространство Рис. 4. Дебит эксплуатационных скважин 


Выводы. Синтезированные нелинейные законы управления м, и›, из, и, 
и; методом АКАР обеспечивают асимптотическую устойчивость движения в 
целом замкнутой системы относительного заданного состояния, определяе- 
мого требованиями процесса управления. Для того чтобы забрать из экс- 
плуатационных скважин 1, 2, 3, 4, 5 м? /сут необходимо понизить уровень 
водоносного горизонта в эксплуатационных скважинах на 0.1-1,0 м. На 
рис.1 видно, что стабилизация переходного процесса (понижение уровня в 
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контрольных скважинах) произойдет примерно через 9 суток после начала 
действия возмущающего воздействия. 


Библиографический список 

1. Атрощенко О.И. Основные положения синергетического подхода. 
Постановка нелинейной проблемы аналитического конструирования аг- 
регированных регуляторов / О.И. Атрощенко // Сб. науч. тр. Сев.-Кав. гос. 
техн.ун-та. - 2008. — № 4. - С.16-19. 

2. Атрощенко О.И. Синергетический синтез системы управления де- 
битом водозаборных скважин методом АКАР / О.И. Атрощенко // Вестник 
Сев.-Кав. гос. техн. ун-та. - 2008. — № 2(15). - С.53-57. 

3. Малков А.В. Синтез распределенных регуляторов для систем 
управления гидролитосферными процессами / А.В. Малков, И.М. Першин. 
— М.: Научный мир, 2007. 

4. Першин И.М. Анализ и синтез систем с распределенными пара- 
метрами / И.М. Першин. — Пятигорск: Рекламно-информационное агентство 
на КМВ, 2007. 

5. Современная прикладная теория управления: Синергетический 
подход в теории управления; под ред. А.А. Колесникова. — Ч.П. —Таганрог: 
Изд-во ТРТУ, 2000. 


Материал поступил в редакцию 4.07.08. 
О.Т. АТКОЗНСНЕМКО 


ЗТМЕКСЕТТСНЕ$КУ $УМТНЕ$1$ ОЕ МАМАСЕМЕМТ$ 
РОК ТНЕ МОМЕТМЕАК ЗРАТТАН-У-МОТ ТМУАКТАМТ ОВЗЕСТ 
ОЕ МАМАСЕМЕМТ 


АррИсаНоп оЁ теод АКАР («Апа4са! деядпттод ог {Пе аддгедае4 геду!аког$ 
оп {Ле 5еЁ пуаНапЕ уаце{ез») этегде{с$ Неопез оЁГ тападетепе Гог зуп{Пе$5 
ОГ |ам/5 о тападетепЕ Бу зраЧаЙу-пепуанаппут поп!пеаг 4$ 1Ьщед оБбесе 
ОГ тападетепе 15 сопу@еге4. Тпе обесЕ ог тападетепе сПоо5ез маег Поп- 
201% {Пе Сепба! °Ке оЁ Не К$оуод$К деро$Ё оЁ тега! мгакегс. 


АТРОЩЕНКО Олег Игоревич (р.1983), аспирант кафедры «Управление и 
информатика в технических системах» Пятигорского государственного тех- 
нологического университета. В 2005 г. окончил Кисловодский гуманитарно- 
технический институт 

Область научных интересов: управление, математическое моделирование, 
автоматика, системы управления с распределенными параметрами. 

Автор шести научных статей. 


251 


